Shear viscosity coefRcient and relaxation time of causal dissipative hydrodynamics in 

QCD 



T. Koidc 

FIAS, Johann Wolfgang Goethe-Universitydt, Ruth-Moufang Str. 



1, 6O438, Frankfurt am Main, Germany 



E. Nakano 

Extreme Matter Institute, GSI, Planckstr. 1, D-64291 Darmstadt, Germany 



T. Kodama 

Institute de Fisica, Universidade Federal do Rio de Janeiro, C. 



P. 68528, 21945-970, Rio de Janeiro, Brasil 



O 
O 

(N 

1—5 
O 

m 



Oh- 

CD 

> 
O 
rn 

O 

o 



X 



The shear viscosity coefRcient and the corresponding relaxation time for causal dissipative hydro- 
dynamics are calculated based on the microscopic formula proposed in [T. Koide and T. Kodama, 
Phys. Rev. E 78, 051107 (2008)]. Here, the exact formula is transformed into a more compact 
form and applied it to evaluate these transport coefficients in the chiral perturbation theory and 
perturbative QCD. It is shown that in the leading order calculation, the causal shear viscosity co- 
efficient rj reduces to that of the ordinary Green-Kubo-Nakano formula, and the relaxation time 
is related to rj and pressure P by a simple relationship, Tt, — rj/P. 



The study of dissipative processes in relativistic heavy 
ion collisions is now one of the important topics to clarify 
quantitatively how closely the matter created there be- 
haves as an ideal fluid. The proper concept of relativistic 
dissipative hydrodynamics is, however, not trivial at all. 
It is by now known that a naive covariant extension of 
the Navier-Stokes equation leads to the problem of rela- 
tivistic acausality and instabilities [l[ . 

An essential key to this question is the presence of 
memory effect characterized by a finite relaxation time 
in the definition of irreversible currents Conse- 
quently, in a relativistic regime, any fluid becomes non- 
Newtonian, that is, the irreversible current is not simply 
proportional to thermodynamic forces. We refer to those 
theories which incorporate this effect as the causal dissi- 
pative hydrodynamics (CDH). 

Similarly to equation of states, transport coefficients 
reflect the properties of the matter. These are inputs 
for hydrodynamic modelings and should be determined 
from a microscopic theory. So far, there are mainly two 
approaches to estimate these coefficients of CDH. One is 
the kinetic approach based on the Boltzmann equation 
[1, 3] ■ The other is the calculation from a duality assum- 
ing the AdS/CFT (Anti-de Sitter space/ conformal field 
theory) correspondence d, However, the applicabil- 
ity of these approaches is not obvious for the physics of 
heavy- ion collisions. The kinetic approach is only appli- 
cable for rarefied gases where the Boltzmann- Gr ad limit 
is satisfied. The AdS/CFT approach predicts the behav- 
ior of a matter described by a conformal field theory in 
its the strong coupling limit, but QCD is obviously not 
a conformal theory. 

The Green-Kubo-Nakano (GKN) formula is another 
possibility which is used to evaluate these transport co- 
efficients in terms of microscopic correlation functions. 
However, for our purpose, there exist the following prob- 



lems. One is that, in this approach, the relaxation time 
has not been expressed with correlation functions. Fur- 
thermore, it is not obvious whether the GKN formula 
for the shear viscosity coefficient is applicable when the 
system exhibits a non-Newtonian nature of irreversible 
currents. Usually Newtonian nature is assumed in the 
derivation of the GKN formula but, as mentioned above, 
this is exactly the key question for the case of a rel- 
ativistic fluid 0, Si- More precisely, the shear viscos- 
ity coefficient appears as a phenomenological parameter, 
that is, the diffusion constant for transverse momentum 
via a ratio to entropy density s in the phenomenological 
Langevin approach. However, this proportionality to the 
entropy density is not always true and depends on its 
definition from a microscopic theory. 

Another important point is that the hydrodynamic 
shear viscosity, for example, appears in a thermal re- 
laxation process which is not a response of the system 
to an external force, but is rather a response to the in- 
homogeneity of the velocity field, that is, the change in 
boundary conditions. Thus a direct application of the 
traditional argument of the linear response theory for the 
calculation of the shear viscosity is not straightforward. 
As a matter of fact, the GKN formula of the shear viscos- 
ity is derived by using, for example, the nonequilibrium 
statistical operator method 

Recently, the microscopic formulae to calculate the 
transport coefficients of CDH was derived by use of the 
projection operator method 0, Bll^l- Although this ap- 
proach is shown to be very powerful, the results obtained 
so far are yet too formal for the practical applications. 
The aim of the present work is two fold: one is to reduce 
the formula obtained in 0, Q into a more compact form, 
and the other is to apply it to the hadronic matter and 
quark-gluon plasma (QGP) to investigate the tempera- 
ture dependence of the causal shear viscosity coefficient 
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r] and the relaxation time r^. 

So far, Ttt has not been discussed enough compared 
to r], but some crucial questions depend on the precise 
behavior of the relaxation time. For example, the behav- 
ior of is directly related to the stability of relativistic 
fluids as discussed in Another question is the value 
of 77/s obtained form the analysis of experimental data 
seems very small, as close, if not smaller than, to the 
lower bound I/Att predicted by the AdS / CFT correspon- 
dence. However, even if 77/ s is not so small, the fluid 
can behave as an ideal fluid when r^r is sufficiently large. 
Thus to know the precise value of r^r is fundamental. 

For simplicity, let us consider a system with shear flow, 
where the fluid velocity points to the x direction with 
finite velocity gradient in the y direction. We further 
assume that the evolution equation of the shear viscos- 
ity is approximately given by the evolution of the two 
variables, T^^ and T°^. Correspondingly, we will em- 
ploy the formula of = 2 form given in Sec. VIII of 
where A^ is the dimension of the projected space. The 
exact time evolution of spatial, off-diagonal components 
of energy- momentum tensor T^" (k, t) in the momentum 
representation (Fourier transform) is formally expressed 



dtTy%ky,t) 

- / dTE22{ky,T)T 
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where Rj,^ = {Ty-{ky),Ty-{~ky)) / {T^-{ky),T^-{-ky)), 
and the inner product {X, Y) represents Kubo's canon- 
ical correlation {X,Y) = d\(i-^TT:[peqe^" Xe~^"Y]. 
Here, H is the Hamiltonian and (3 is the inverse of tem- 
perature \/T. The memory function S22(fcy,t) is defined 
in Eq. (78) of j!]. The last term of the Eq. ^ is often 
called the noise term, and related to the memory term 
through the fiuctuation-dissipation theorem of the second 
kind (3,il. 

To break the time-reversal symmetry, we introduce the 
time-convolutionless (TCL) approximation0,@|, 

dtTy^{ky,t) = -iky{e + p)Ruy{ky,t) - —Ty^{ky,t), 

Tky 

(2) 

where we have introduced a function related to the relax- 
ation time TjT ^ ~ /q°° drS22 {ky,T), and e , P and u are the 
energy density, pressure and fluid velocity, respectively. 
The noise term is not necessary for the present discus- 
sion and dropped out. The TCL approximation is equiv- 
alent to assume the exponential ansatz for the memory 
of Ty^{ky,t). It should be emphasized that, even after 
the TCL approximation, the hysteresis of T^^ is still re- 
served in its time evolution equation and this corresponds 
exactly to the memory effect in the usual CDH[2[. In 
this derivation, we have used the following replacement 
j^Ox — ^ P)u^{ky,t), which is justified near the local 
rest frame. 



On the other hand, the linearized phe- 
nomenological equation for the shear viscos- 
ity defined by projection to the traceless part, 
T^t^-^ = i (A^"A'''3 + A'^'^A''" - 2^A'>'A"'3) To,p with 
A^'^ = g^'^ — u^u'^ , is given by 



r.^^^'^'iky) + ^"''{ky) = -^[zky)u^{ky) (3) 

at the local rest frame. Comparing the above equation 
with Eq. (I2|), we identify the rj and Q: 

= i.^™n '^^y' ^ = J"^^n(^ + P)RkyTky ■ (4) 

The memory function can be cast into a simple form 
in the low momentum limit: 



(5) 



where X22{ky,s) is the Laplace transform with respect 
to t of the following correlation function. 



X22iky,t) 



iTy-iky,t),Ty%-ky,o)) 

{TV^{ky,0),Ty^{-ky,0))' 



(6) 



Substituting this equation into the definitions ([3]), we ob- 
tain = Xi'2{0, 0) and r] = {e + P)RoX^2iO, 0). 

For later convenience, we re-express the correlation 
function A'^2(0,0) with the definition of the GKN shear 
viscosity coefficient Q as 



where 
Vgkn 



10 



Vgkn 



d3x(7r^^(x,0),7r^^(0,0)) 



(7) 



dti 



dr(^"''(x,t)^„^(0,r)),e 
(a,/? = 0,1, 2, 3). (8) 



Here the suffix ret denotes the retarded Green function, 
and there is difference by factor 2 from Q because of the 
difference of the definition. 
Finally, -q and are given by 



1] 



Vgkn 



P{e + P) /32/d3x(j^(x,0), J-(0,0))' 

TV _ TlGKN 

13 ~ f3^J d3x(7r^y(x,0),7r^!'(0,0))' 



(9) 
(10) 



where we have introduced the energy current, J^(x, t) = 
u^T''''{yi,t) which is T°''(x,t) at the local rest frame. 
Note that ri/{f3{e + P)) is nothing but rj/ s dX vanishing 
chemical potential. These expressions are still most gen- 
eral ones within a framework where the hydrodynamic 
description is meaningful and are the first main results 
of the present work. 
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Concerning the causality problem in CDH, we intro- 
duce a ratio r\l(jT^{e + P)), which is, in our case, given 
by i?o itself. Note that the propagation speed of signal 
should be slower than the speed of Hght. For the 3+1 di- 
mensional viscous fluid to be causal in the shear channel, 
the following relation should be hold [H, 



1] 3 , 
< -fl 



(11) 



where is the ordinary sound velocity = dP/de. In 
particular, is 1/3 in the massless ideal gas case and 
this r] — ration should be smaller than 1/2. 

Now we apply our formula to QCD related problems. 
In the hadronic phase at low baryon density and high 
temperature, the dominant behavior comes from the pi- 
ous. The interaction of pions can well be described by 
the chiral perturbation theory which is a low energy ef- 
fective theory of QCD. To estimate the GKN shear vis- 
cosity coefficient riQ^N, we adapt the leading order result 
given in fllj . At the high temperature regime far above 
the deconfinement scale, we assume that the dominant 
behavior of the QGP is given by gluon, and use rjcKN 
obtained to the leading order in the perturbative QCD 
(pQCD) calculation for iV/ = in [l2|. We take the run- 
ning coupling constant to have the same parametrization 
as |13| . 

In order to estimate i] and r^r to the leading order in 
couplings, it is enough to calculate the fluctuations of J'' 
and TT^" to the lowest order. Then we obtain 



j d=x(.7'(x), J'(0)) = ^P'",(l + 

(12) 

+ /3(1 + rip)^ Up 
(13) 



d3x(7r^y(x),7r^«(0)) 



g I d^p (p2)2 f 1 



15 7 (27r)3 /3P2 
_ P 



where Up = (e^^^ — 1)^^ is the Bose-Einstein distribution 
function. The statistical factor 5 is 3 for pions and 16 for 
gluons. In the derivation of Eq. ([13]), we have subtracted 
the vacuum contribution. Here e and P are the energy 
density and pressure of the ideal Bose gas, respectively. 

By substituting them into the definitions Eqs. ^ and 
(fTU)) . we obtain the leading order results. In this approx- 
imation, we found that the causal shear viscosity coeffi- 
cient is given by rj = rjcKNi and the relaxation time (jlO[) 
reduces to 



Ttt = tigkn/P = 11 /P- 



(14) 



See [111, [ij, llJ] for the behavior of tjgkn- 

We show the temperature dependence of to /? ratio 
in Fig. [T] The behavior is very similar to that of 77/s: a 
decreasing function of temperature in the hadron phase 




FIG. 1: The temperature dependence of the relaxation time 
to (3 ratio with /3 = 1/T. The solid and dashed lines rep- 
resent our result and that from AdS/CFT correspondence, 
respectively. 
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FIG. 2: The temperature dependence of TTrP/rj. The solid, 
dashed and dot-dashed lines correspond to the result of our 
formula, the moment method and AdS/CFT correspondence, 
respectively. 



and an increasing one in the QGP phase, exhibiting min- 
imum near the phase transition temperature Tc- The 
dashed line denotes the prediction from the AdS/CFT 
correspondence t^/P = (2 - ln2)/(27r) for TV = 4 SYM 
0. However, because of the weak tempcratrure depen- 
dence of T,r//3 in the QGP phase, Ttt itself is a decreas- 
ing function for both of the hadron and QGP phases and 
show discontinuity near Tc, that is, becomes minimum 
(maximum) in the hadron (QGP) side near Tc- Thus the 
behavior of the shear viscosity is insensitive for the rapid 
change of the fluid velocity above Tc because of the large 
Ttt, meanwhile it is more sensitive below Tc. In this sense, 
the matter created in relativistic heavy-ion collisions can 
be close to the ideal fluid in different two ways, one is 
because of the small ry and r^r which may be realized in 
the hadron phase and the other is the small rj but large 
in the QGP phase. 

The relaxation time is calculated also from the 
relativistic Boltzmann equation using Grad's moment 
method with the 14 moments approximation [i, 15|. In 
order to compare with our result, the temperature de- 
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FIG. 3: The temperature dependence of the ratio. 



pendence of Ttt to rj ratio in the unit of pressure P is 
shown in Fig. [2l The sohd hne denotes the behavior 
of our result, which should be one in our leading or- 
der calculation. The dot-dashed line denotes the result 
from the AdS/CFT correspondence, (2-ln2)/2 ^. The 
dashed line represents the result obtained from the mo- 
ment method for the pion (Bose-Einstein) gas without 
phase transition This figure shows that our theory 
gives a different result from the "14" moments calcula- 
tion. In the Navier-Stokes limit (Newtonian fluid), it is 
known that the momentum method with the " 13" mo- 
ments approximation and the famous Chapman-Enskog 
procedure are consistent for the shear viscosity. How- 
ever, Ttt requires expansions in higher moments and, to 
see the relation of our formula and the moment method, a 
more careful comparison should be done. It is also worth 
mentioning that there is a attempt to calculate these co- 
efficients without using the moment method [lij . 

As the AdS/CFT correspondence predicts that the 
minimum of tjgkn/s has a lower bound l/(47r) our 
rj/s and r^r would have a lower bound, unless the fluc- 
tuation of tt'^'^ diverges. As a matter of fact, to be con- 
sistent with the causality condition (jlip . is somehow 
correlated with 77, and t^//3 cannot be much smaller than 
r]/{f3{e + P)). In Fig. [3] ,we show the temperature depen- 
dence of the 77- ratio. In our leading order calculaiton, 

and shows a 



this ratio is nothing but ^ {!^+p) ~ (e+p) ' 
non-trivial temperature dependence only at lower tem- 
perature than the pion mass and finally converges to the 
massless ideal gas limit 1/4 where e = 3P. One can eas- 
ily see that the result satisfies the relativistic causality 
condition (fTT|) . It should be emphasized that relativistic 
fiuids become unstable if the causality condition is not 
satisfied [l[. To see the consistency of CDH, it is neces- 
sary to investigate this ratio more carefully. 

In summary, we proposed the compact definitions of 
the transport coefficients of CDH, and calculated them 
in the chiral perturbation theory and pQCD. We found 
that, in the leading order calculation, the causal shear 
viscosity coefficient 77 is reduced to that of the GKN shear 
viscosity coefficient rjcKN- Although intuitive, this result 



is not trivial a priori, because the irreversible current is 
not simply proportional to the thermodynamic force in 
CDH. In fact, r^r does not vanish in our calculation, and 
there is a simple relationship, r^r = rj/P. This relation is 
not the same as that obtained in the 14 moment calcu- 
lation, and the physical origin of this discrepancy should 
be clarified. So far, the improvements of the moment 
method has been discussed in diverse ways [l7| . Our re- 
sult may give a milestone in such a development. 

It is thus very interesting to ask how these are modified 
when the calculation is implemented beyond the leading 
order. The simple relationship obtained here would not 
be satisfied in the strong coupling limit. In fact, from 
AdS/CFT results, we can show tj/t^ = rs/(4 - 2 In 2), 
which is larger than the equation of state of the mass- 
less ideal fluid, P = Ts/A, that is, P ^ Tj/r.,^ as was 
shown in Fig. [2l Thus we expect that the simple rela- 
tion 77 = rjGKN and/or = rj/P calculated from the 
exact formulae ([9|) and (jlOp would acquire a non-trivial 
temperature dependence once going beyond the leading 
order result and in the strong coupling limit. 

T. Koide thanks to X. Huang and S. Pu for useful 
comments. This work was (financially) supported by 
FAPERJ, CNPq and the Helmholtz International Cen- 
ter for FAIR within the framework of the LOEWE pro- 
gram (Landesoffensive zur Entwicklung Wissenschaftlich- 
Okonomischer Exzellenz) launched by the State of Hes- 
sen. 



[1] 
[2] 

[3] 



[4] 

[5] 
[6] 



[7] 
[8] 
[9] 
[10] 

[11] 

[12] 
[13] 
[14] 
[15] 
[16] 
[17] 



G. S. Denicol et al, J. Phys. G 35, 115102 (2008). 
See, for example, G. S. Denicol et al, J. Phys. G 36, 
035103 (2008), and references therein. 
See, for example, R. Baier, et al., Phys. Rev. C 73, 
064903 (2006); A. Muronga, Phys. Rev. C76, 014910 
(2007); Z. Xu and C. Greine r, Phys. Rev. Le tt. 100, 
172301 (2008); B. Betz et al, larXiv:0812. 14401 [nucl-th] 
and references therein. 

W. Israel and J. M. Stewart, Ann. Phys. (N.Y.) 118, 341 
(1979). 

P. Kovtun et al, Phys. Rev. Lett. 94, 111601 (2005). 
R. Baier et al, JHEP 0804 (2008) 100; M. Natsuume and 
T. Okamura, Phys. Rev. D77, 066014 (2008); Erratum 
ihid. D78, 089902 (2008). 
T. Koide, Phys. Rev. E75, 060103(R) (2007). 
T. Koide and T. Kodama, Phy. Rev. E 78, 051107 (2008). 
A. Hosoya et al, Ann. Phys. (N.Y.) 154, 229 (1984). 
R. Zwanzig, Noneqmhbrium Statistical Mechanics (Ox- 
ford Univerasity, New York, 2004). 

J.-W. Chen and E. Nakano, Phys. Lett. B647, 371 
(2007). 

P. Arnold, et al., JHEP 05, 051 (2003). 

L. P. Cernai et al, Phys. Rev. Lett. 97, 152303 (2006). 

R. Lacey et al., Phys. Rev. Lett. 98, 092301 (2007). 

M. Prakash, et al., Phys. Rep. 227, 321 (199 3). 

M. A. York and G. D. Moore, [arXiv:0811.0729 [hep-ph]. 

See, for example, A. G. Bezerra et al. Continuum Mech. 



Thermodyn. 6, 149 (1994). 



